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J> \ Results from the Wroclaw Monte Carlo neutrino generator of events 

are reported. Predictions for charged hadron multiplicities, neutral pion 
and strange particle production are presented and compared with available 
\ data. 
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\D ' PACS numbers: 13.15.+g, 17.87.Fh, 14.20.Jn 
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Q H ' 1. Introduction 

Ph! The aim of this paper is to present some details about the performance 

of the neutrino Monte Carlo generator of events developed at the Wroclaw 
University. It differs from other generators (e.g. NEUGEN, NUANCE, 
NEUT) PP in the treatment of the resonance region. Our generator con- 
^ ■ tains an explicit A resonance excitation model but contributions from more 

massive resonances are absent. It is assumed that an average description 
of the cross section coming from those resonances is sufficient. In fact, in 
the low energy BNL data only the A peak is clearly seen in all single pion 
production (SPP) channels on free nucleon targets [2]. In neutrino inter- 
actions with nucleus targets the Fermi motion is supposed to average the 
cross section contributions from other then the A resonance peaks. 

The performance of the Wroclaw generator in the SPP channels has been 
discussed elsewhere (3]. The agreement with available data is satisfactory. 
In order to cover the entire allowed kinematical region it is necessary to 
apply the DIS formalism to produce the events also for values of invariant 
hadronic mass W < 2 GeV. Two problems arise then: (i) how to model the 
structure functions in order to reproduce the correct value of the inclusive 
cross section and (ii) how to produce the final states. 



* Presented by J. A. Nowak at 20th Max Born Symposium Nuclear effects in the neu- 
trino interactions, December 7-10, 2005, Wroclaw, Poland 



(1) 



2 




10 100 10 100 

W 2 (GeV 2 ) W 2 (GeV 2 ) 

Fig. 1. Average charged hadron multiplicity in neutrino scattering off proton and 
neutron targets. Data points are taken from |12l ITU] . 

The generally accepted way to describe the structure functions is to use 
Bodek-Yang low Q 2 modifications, which have been introduced on the basis 
of the electron scattering data [I]. Recently the neutrino scattering data 
from the CHORUS and NuTeV experiments on nuclear targets (iron and 
lead) have become available 5 and further progress and cross checks will 
be possible. In order to produce final hadronic states several strategies can 
be adopted. In our approach we assume that interactions take place on 
separate constituent of the nucleon. The first step of fragmentation is per- 
formed by our generator and of remaining quark-diquark system by means 
of PYTHIA6 routines based on the LUND model [B]. Several parameters 
of PYTHIA6 were fine-tuned to get a good agreement with data [7]. In 
effect it is possible to produce DIS events for small values of the invari- 
ant hadronic mass down to the threshold for the single pion production 
W > Wthr = M + m n . The low W DIS events are used in order to model 
the non-resonant background [S]. 

The description of the whole generator is presented in j^J. In this paper 
we will present some details about the performance of generator's hadroniza- 
tion routines. 

2. Results 

Average multiplicities (n c h) of charged hadrons have been measured in 
several experiments and the typical dependence 

(n ch ) = a + b -In W 2 / GeV 2 , (1) 

has been observed. In ref. ^U] (it contains a missprint in the value of a in 
P|) ) the following fits have been found: 
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Fig. 3. Average charged hadron multiplicity in antineutrino scattering on proton 
and neutron. Data points are taken from J2| an d the fit was found in 
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a = 0.50 ±0.08 b =1.42 ±0.03 for CC vp -» fiX ++ , (2) 
a = -0.20 ±0.07 b = 1.42 ±0.03 for CC vn -> . (3) 

In the experiment ^Uj results for scattering off proton and neutron target 
were extracted from the Fermilab 15-foot deuterium bubble chamber data. 
Average neutrino energy was about 50 GeV and multiplicity distributions 
were studies in the invariant mass range 1 < W < 15 GeV. 

In the paper quite different values of the parameter a nd b were 
found: 

a = -0.05 ±0.11 b =1.43 ±0.04 for CC vp -»• fj,X ++ . (4) 

In the experiment ^2] multiplicities of hadrons were obtained from the 
hydrogen bubble chamber BEBC and the investigated range of the invariant 
hadronic mass was 3.5 < W < 10 GeV. 

The results from the Wroclaw generator are shown in Fig. ^ It is 
seen that the average multiplicities are much lower then the data from ^UJ 
while agreement with data from |12j is very good. The difference between 
two measurements is probably due to rescattering inside deuterium [TP] . 
We verified that very similar to ours results are produced by the original 
PYTHIA6 code. 

The distribution of multiplicities agrees well with the KNO model and 
is independent of the value of W. It means that the probability to observe 
n charged hadrons is 



f(n, (n ch )) = F , 



(5) 
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Fig. 6. Average p° production multiplicities in (anti-) neutrino scattering off proton. 
Data points are taken from |2J 



where F is the universal scaling function. 

In Fig. [2 the distributions of multiplicities obtained from the Wroclaw 
generator are shown and compared with the data from ^Uj- We see that 
the pattern of distributions predicted by the MC generator agrees with the 
data but higher multiplicities are systematically underestimated as could be 
expected from the results for (n c ^). 

In Fig. |21 we show the average charged hadron multiplicities for antineu- 
trino reactions. The fits obtained in are 

a = 0.02 ± 0.20 b = 1.28 ± 0.08 for CC up -> fi + X° (6) 
a = 0.80 ±0.09 b = 0.95 ±0.04 for CC un^fi + X~. (7) 

The agreement of generator's predictions with the data is satisfactory. 

It is also important to investigate NC reactions. Here the data is avail- 
able with the average charged hadron multiplicities [TB] : 

a = 1.61 ±0.16 b = 0.99 ±0.05 forNC vp -> uX + (8) 
a = 1.65 ±0.20 b = 0.95 ±0.07 for NC vp^>vX + . (9) 

The comparison of MC predictions with the data is shown in Fig. |3] and 
the clear disagreement is seen. We were not able to obtain the original 
paper and data. The presented fits are taken from [T2|. As far as we know 
the NC data was never published as an official WA21 collaboration's paper. 

We did not find the experimental data for charged pion multiplicities. 
One can use the existing data from anti-muon scattering and estimate the 
charged pion multiplicities in neutrino reactions but such reconstruction 
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Fig. 7. Cross sections for channel vn — > /.i~K + A. Data points are taken from 

QUE] 

carries especially for lower values of W an uncertainty which is difficult to 
evaluate. 

The precise data exists for neutral pion production on the proton target. 
They are again given in the form of a fit to the general formula (QQ): 

a = 0.14 ±0.26 b = 0.50 ±0.08 for CC up vr° production (10) 
a = 0.17 ±0.42 b = 0.53 ±0.13 for CC up vr° production. (11) 

In the Fig. [5] we show our MC predictions for ir° production. The 
agreement with the data is excellent. 

In the Fig. |f)]we show our generator's predictions for the p° production 
in up and up scattering. The data is taken from |13j . In the case of neutrino 
the agreement with data is satisfactory, but for anti-neutrino scattering the 
cross section for p° production is overestimated. 

Another good test of the hadronization routines of the generator is pro- 
vided by predictions for the strangeness production. 

At low anti-neutrino energy strange particles can be produced in quasi- 
elastic processes but we do not discuss them here. In Fig. |7|the generator's 
predictions for the cross section for un — > fi~K + A reaction are shown. The 
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data was taken from |14| 115] and it does not seem to be consistent (notice 
the log scale). The MC results are in general agreement with the data. 

Finally, in the Fig. |H] we show the cross section for neutral strange 
particles K° and A . Data points are taken from [1711181 HH1I2U] . We notice 
the good agreement of the generator's predictions with the data. 

3. Conclusions 

The predictions of the Wroclaw MC generator of events are in satisfac- 
tory agreement with the data. Nevertheless additional effort is necessary to 
fine tune some of the free parameters of the generator in order to improve 
it's performance. The work on the nuclear effects module of the generator 
is in progress. 

Acknowledgements 

The authors would like to thank C. Juszczak for many conversations and 
for participating in this project at various stages. The authors were sup- 
ported by KBN grant 105/E-344/SPB/ICARUS/P-03/DZ211/2003-2005 and 
by Wroclaw University grant 2595/W/IFT. 



REFERENCES 

[1] Proceedings of the series of Nulnt workshops: Nucl. Phys. B (Proc. Suppl.) 

112 (2002); ibid. 139 (2005) 
[2] T. Kitagaki et al, Phys. Rev. D 34 (1986) 2554. 

[3] C. Juszczak, J. A. Nowak and J. T. Sobczyk, Simulations from a new neutrino 
event generator, arXiv:hep-ph/0512365 



8 



[4] A. Bodek and U. K. Yang, Nucl. Phys. Proc. Suppl. 112 (2002) 70. 

[5] |http: / /neutrinooscillation.org/| |http : / / www-e8 1 5 . f nal . gov /| . 

[G] T. Sjostrand, et al., Comput. Phys. Commun. 135(2001) 238 

[7] J. A. Nowak, Neutrino event generator, to be published in Physica Scripta T 

[8] J.T. Sobczyk, J.A. Nowak, K.M. Graczyk, Nucl. Phys. B (Proc. Suppl.) 139 
(2005) 266 

[9] J.A. Nowak Construction of the Monte Carlo generator of neu- 
trino interaction, PhD Thesis (in Polish), will be available at 
http: / /wng. ift.uni.wroc.pl/theses. html 

[10] D. Ziemihska et. al., Phys. Rev. D27 (1983) 47 

[11] S. Barlag et al., Z. Phys. Cll (1982) 283 

[12] H. Grassier et al., Nucl. Phys. B223 (1983) 269 

[13] H. Grassier et al., Nucl. Phys. B 272 (1986) 221 

[14] N. J. Baker et al, Phys. Rev. D 24 (1981) 2779. 

[15] D. Son et al, Phys. Rev. D 28 (1983) 2129. 

[16] E. Hoffman, Multiplicity in NC vp Dp scattering, PhD Thesis cited in 
N. Schmitz, Hadron Production in High-Energy Neutrino and Anti-Neutrino 
Collisions, preprint MPI-PAE/Exp-El-193. 

[17] N.J Baker et al. Phys. Rev D34 (1986) 1251 

[18] P. Bosetti et al. Nucl. Phys. B209 (1982) 29 

[19] D. DeProspo et al., Phys.Rev. D50 (1994) 6691 

[20] G.T. Jones et al. Z Phys. C 57 (1993) 197 

[21] P. Allen et al. Nucl. Phys. B 194 (1982)373 



